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Dynamics of Particle Production in Relativistic Nuclear Collisions 
K. Tuominen a 

a NORDITA, Blegdamsvej 17, 2100 Copenhagen 0, Denmark 

Saturation models for particle production in relativistic nuclear collisions are discussed. 
In particular, I show that the predictions from the high density QCD for the qualitative 
shape of dN/dy are very sensitive to the form of the unintegrated gluon distribution. 

1. Introduction 

; After first years of running, RHIC has provided a lot of interesting data [1]. The very 
basic observable, the number of particles at central unit of rapidity, seems to indicate that 
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saturation models describe well the center of mass energy dependence as well as the cen- 



trality dependence. The models I consider in this talk are: The pQCD+saturation -model 



[2] and the high density QCD calculation by Kharzeev and Levin (KL) [3]. Although both 
models agree with data and each other at central rapidity, they show qualitative differ- 



ences away from midrapidity. We will trace the origin of this discrepancy and see how it 

could be improved upon and, along the way, we shall briefly remark how these two models 

D ! would predict the second important observable, the transverse energy. 
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2. Models and results 

In the pQCD+saturation model the multiplicity of initially produced gluons is evalu- 
ated under the assumption of collinear factorization and including all quanta above the 
saturation scale p s which is obtained as a self consistent solution of the saturation condi- 
tion 
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After solving for p s one obtains the scaling laws [2] 
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for the initial gluon multiplicity and transverse energy at y = and, assuming ideal ID 
expansion, for the multiplicity and transverse energy of hadrons at y — 0. In (2) T c is 180 
MeV and evolution in the hadronic phase has been neglected as this is compensated for 
by the development of the flow. 
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The KL calculation also predicts powerlike growth of the number of particles in central 
rapidity with y/s [3]: 



N ~ N 



part I 



A 2 



(3) 



where A = 0.25 is related to the small Bjorken-a; growth of the gluon structure function. 
The power is slightly smaller than the one in the pQCD+saturation calculation. The figure 
1 shows the results from these two models and one sees that in the presently available 
energy range the models are indistinguishable at central (pseudo)rapidity. Extrapolation 
to LHC energies leads to a wider range of predictions and not all models are distinguishable 
even there. For the pQCD+saturation curve an effective value of 178 for A was used, as 
this corresponds to the 6% centrality cut of the data [4]. A crude estimate of E± in the 
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Figure 1. Multiplicity in the central unit of Figure 2. dN/dy from the two saturation 



pseudorapidity. Data is from PHOBOS and 
dashed lines show the log(s) (straight) and 
log 2 (s) (curved) growths. 



models. Normalization of the KL curve is 
such that after transforming, with an over- 
all Jacobian factor, from y to 77 it agrees 
with data. Figure is from [5]. 



pQCD+saturation model can be obtained from Eq.(2): at y/s = 130 GeV E± = 520 GeV 
and at y/s = 200 GeV E± = 620 GeV. Transverse expansion effects increase the above 
numbers a little [4]. PHENIX data at y/s = 130 is E ± = 578 GeV [6]. 

The KL calculation has been shown to reproduce pseudorapidity distributions of hadrons 
at y/s = 130 GeV well [3]. From the theoretical point of view it is more instructive to look 
at the rapidity distributions of initial gluons, see Fig. 2. The pQCD+saturation model 
leads to a broad gaussian which, after transforming to pseudorapidity and comparing with 
data, overshoots at large rapidities. If saturation dynamics lead to flat dN/dy near y ~ 0, 
this behaviour has to change around y ~ 2/beam/2 and the saturation dynamics must be 
replaced by some other fragmentation region dynamics. The KL result is qualitatively 
very different: it has a discontinuity in the first derivative of dN/dy at y = and is 
exponentially suppressed away from it. 
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2.1. High Density QCD methods 

The result of KL is based on the GLR equation [7] : 

= — / dp\^r [ P± dklfafa, k 2 ± )(j) A (x 2 , (p - k)\), (4) 
which requires an ansatz for (f) A (x, kj_). The ansatz used by KL is effectively 



M*, k\) = ^- 2 ^[9(Ql A - kl) + e^6(ki - Q sA )\, (5) 



where parametrically e ~ 0(a 2 ) and sets the relative normalization of the saturated and 
perturbative parts of the unintegrated gluon distribution. KL choose e = and regard 
the tail as a small correction. With ansatz (5) one obtains to leading logarithmic accuracy 

Setting e = the KL result is reproduced. With e = 1 and A = 0.25 one sees that the 
term 2e[. . .] originating entirely from the tail of the distribution and the first term from 
saturation region are roughly equal. For the transverse energy the tail is more important: 

Of course it is plausible that the relative normalization is not given by e = 1 but by some 
smaller value, since e ~ 0(a 2 ), which would validate the exclusion of the tail. 

2.2. Reshaping dN/dy 

The ansatz (5) is probably too simple and one should try for example the one from [8]: 

4> A {x, k\) ~ / —^\\ - e -kQU*)* 2 H4/z 2 +V). ( 8 ) 
J az 2 

Using this in (4) one finds that the form of dN/ dy is 

^ - 7 [6 + 1 (2Ay) 2 - — (2A# + ...]. (9) 

dy a(Q s>A ) sinh(Ay) 1 18 V y; 1800 V yj 1 K ' 

This is very different from the one obtained with ansatz (5) and is closer to a broad 
gaussian as can be seen from the figure 3, which shows all of the discussed multiplicity 
distributions. 

One needs also to take into account the large x behaviour of the gluon distribution, 
xG ~ (1 — re) 4 , not contained in Eq.(8); whether this is sufficient to make agreement with 
data remains unclear at the moment. However near y ~ 0, one expects result (9) to be 
dominated by small x part of <f) A , and one might try to transform to pseudorapidity and 
compare with data. For KL this is done by an overall factor and for pQCD+saturation 
-model by assuming exponential p_i_-spectra for pions. For details, see [3,5]. From Fig. 4 
one sees that the transformation of (9) with the overall factor leads to a very large dip 
around the central pseudorapidity. Hence, one should rather use the p^-distributions to 
carry out the transformation. Probably the effects of the large x behaviour of the gluon 
distribution should be included already at y ~ 0, too. 
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Figure 3. dN/dy from the different calcula- Figure 4. dN/di] from the different calcu- 
tions discussed in the text. lations discussed in the text. 

3. Conclusions 

Particle production over the whole experimentally accessible rapidity range has been 
investigated using saturation models. While these models lead to very similar results 
at central unit of rapidity, they seem to differ at nonzero rapidities. This difference was 
suggested to originate from the choice for the unintegrated gluon structure function in the 
KL calculation, and a different ansatz was shown to lead to a gaussianlike distribution 
with a width comparable to the pQCD+saturation -model result. The quality of the 
approximations such as the neglect of the tail of the distributions and transformation 
from y to rj with an overall Jacobian factor, was shown to be model dependent. 
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